Classical capacity of bosonic broadcast communication 
and a new minimum output entropy conjecture 
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Previous work on the classical information capacities of bosonic channels has established the 
capacity of the single-user pure-loss channel, bounded the capacity of the single-user thermal-noise 
channel, and bounded the capacity region of the multiple-access channel. The latter is a multi-user 
scenario in which several transmitters seek to simultaneously and independently communicate to a 
single receiver. We study the capacity region of the bosonic broadcast channel, in which a single 
transmitter seeks to simultaneously and independently communicate to two different receivers, ft 
is known that the tightest available lower bound on the capacity of the single-user thermal-noise 
channel is that channel's capacity if, as conjectured, the minimum von Neumann entropy at the 
output of a bosonic channel with additive thermal noise occurs for coherent-state inputs. Evidence 
in support of this minimum output entropy conjecture has been accumulated, but a rigorous proof 
has not been obtained. In this paper, we propose a new minimum output entropy conjecture that, if 
proved to be correct, will establish that the capacity region of the bosonic broadcast channel equals 
the inner bound achieved using a coherent-state encoding and optimum detection. We provide some 
evidence that supports this new conjecture, but again a full proof is not available. 

PACS numbers: 03.67.Hk, 89.70.+C, 42.79.Sz 
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I. INTRODUCTION 

The past decade has seen several advances in evalu- 
ating classical information capacities of several impor- 
tant quantum communication channels Despite 
these advances | lj . exact capacity results are not known 
for many important and practical quantum communica- 
tion channels. Here we extend the line of research aimed 
at evaluating capacities of bosonic communication chan- 
nels, which began with the capacity derivation for the 
input photon-number constrained lossless bosonic chan- 
nel The capacity of the lossy bosonic channel 
was found in Q, where it was shown that a modulation 
scheme using classical light (coherent states) suffices to 
achieve ultimate communication rates over this channel. 
Subsequent attempts to evaluate the capacity of the lossy 
bosonic channel with additive Gaussian noise [5[ led to 
a crucial conjecture on the minimum output entropy of 
a class of bosonic channels @. Proving that conjecture 
would complete the capacity proof for the bosonic chan- 
nel with additive Gaussian noise, and it would show that 
this channel's capacity is achievable with classical-light 
modulation. More recent work that addressed bosonic 
multiple-access communication channels [7[ revealed that 
modulation of information using non-classical states of 
light is necessary to achieve ultimate single-user rates in 
the multiple-access scenario. 

In the present work, we study the classical informa- 
tion capacity of the bosonic broadcast channel. A broad- 
cast channel is the congregation of communication me- 
dia connecting a single transmitter to two or more re- 
ceivers. In general, the transmitter encodes and sends 
out independent information to each receiver in a way 
that each receiver can reliably decode its respective in- 
formation. We will show that when coherent-state en- 
coding is employed in conjunction with coherent de- 



tection, the bosonic broadcast channel is equivalent to 
a classical degraded Gaussian broadcast channel whose 
capacity region is known, and known to be dual to 
that of the classical Gaussian multiple-access channel 
Q. Thus, under these coding and detection assump- 
tions, the capacity region of the bosonic broadcast chan- 
nel is dual to that of the multiple-access bosonic chan- 
nel with coherent-state encoding and coherent detection. 
To treat more general transmitter and receiver condi- 
tions, we use a limiting argument to apply the degraded 
quantum broadcast-channel coding theorem for finite- 
dimensional state spaces [9] to the infinite-dimensional 
bosonic channel with an average photon-number con- 
straint. We consider the two-receiver case in which Al- 
ice (A) simultaneously transmits to Bob (B), via the 
transmissivity r\ > 1/2 port of a lossless beam splitter, 
and to Charlie (C), via that beam splitter's reflectivity 
1 — 77 < 1/2 port, using arbitrary encoding and optimum 
measurement with an average photon number N at the 
input. Given a new conjecture about the minimum out- 
put entropy of a lossy bosonic channel, we show that the 
ultimate capacity region is achieved by coherent-state en- 
coding, and is given by 

Rb < 9(r)0), Rc < g((l-v)N)-g((l-v)PN), (1) 

for < (3 < 1, where g(x) = (x + 1) log(.T + 1) — a;log(a;) 
is the von Neumann entropy of the Bose-Einstein distri- 
bution with mean x. Interestingly, this capacity region is 
not dual to that of the bosonic multiple-access channel 
with coherent-state encoding and optimum measurement 
that was found in [7|. 

The remainder of this paper is organized as follows. 
Section II gives a brief introduction to the capacity region 
of classical broadcast channels. In Sec. Ill, we describe 
some recent work on the capacity region of the degraded 
quantum broadcast channel Q. In Sec. IV, we introduce 
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the noiseless bosonic broadcast channel model, and de- 
rive its capacity region subject to a new minimum output 
entropy conjecture. In Sec. V we compare the capacity 
region obtained in Sec. IV with the classical Gaussian 
broadcast channel results that apply for coherent-state 
encoding and coherent (homodyne or heterodyne) detec- 
tion. We also show that a recent duality result between 
capacity regions of classical multiple- input, multiple- 
output Gaussian multiple-access and broadcast channels 
Q does not hold for bosonic channels with coherent- 
state encoding. Discussion of bosonic-channcl minimum 
output entropy conjectures, and evidence supporting the 
conjecture associated with the bosonic broadcast chan- 
nel, will be given in Appendix A. 

II. CLASSICAL BROADCAST CHANNEL 



Degraded broadcast channels were first studied by Cover 
[ig |. who conjectured that the capacity region for Alice 
to send independent information to Bob and Charlie at 
rates Rb and Rc respectively over such a channel is the 
convex hull of the closure of all (Rb, Rc) satisfying 

Rb < I(A;B\T) (7) 
Rc < I(T;C), (8) 

for some joint distribution PT(t)PA\T( a I t)PB,c\A(P, 7 I 
a). Here, I(X;Y) denotes the Shannon mutual in- 
formation between ensembles X and Y, and T is 
an auxiliary random variable with cardinality \T\ < 
min{|^4|, \B\, \C\}. The achievability of the above capac- 
ity result was proved by Bergmans [m, and Gallager 
proved the converse fl2l ]. 



A two-user discrete, memoryless broadcast chan- 
nel is modeled by a classical probability distribution, 
Pb,c\a(P,1 I a ), where a, (3, and 7 are drawn from 
Alice's input alphabet A, and Bob and Charlie's out- 
put alphabets B, C respectively. A broadcast chan- 
nel is said to be memoryless if successive uses of the 
channel are independent, i.e., Ps n ,c n \A n [fl n , 7™ | oi n ) — 
Tli=iPB,c\A(/3i, li I o*i) is the transition distribution for 
n channel uses. A ((2 ,iRb , 2 nRc ), n) code for a broadcast 
channel consists of an encoder 

a n . 2 nR B x 2 ni?. c _> 

and two decoders 

W B : B n -> 2 nRs (3) 
W c : C n -> 2 nRc . (4) 

in) 

The probability of error P e is the probability that the 
overall decoded message does not match the transmitted 
message, i.e., 

P {n) = P (w B {B n ) ± W B or W c (C n ) ^ W c ), (5) 

where the messages, (Wb, Wc), that are sent to Bob and 
Charlie, respectively, are assumed to be uniformly dis- 
tributed over the 2 uRb x 2 nRc possibilities. A rate pair 
(Rb, Rc) is said to be achievable, for the broadcast chan- 
nel, if there exists a sequence of ((2 uRb , 2 nRc ), n) codes 

in) 

with Pe — > as n — > 00. The capacity region of the 
broadcast channel is the closure of the set of achievable 
rates. 

Determining the capacity region of a general broadcast 
channel is still an open problem. The capacity region 
is known, however, for degraded broadcast channels, in 
which one receiver (say C) is "downstream" from the hrst 
receiver (say B), so that C always receives a degraded 
version of what B observes. In other words, A — > B — > C 
is a Markov chain, so that there exists a distribution, 
p(j I (3), such that 

Pc\Ah\a)=J2PB\A(P\a)p( 7 \/3). (6) 



III. QUANTUM DEGRADED BROADCAST 
CHANNEL 

A quantum channel Na-b from Alice to Bob is a trace- 
preserving completely positive map that transforms Al- 
ice's single-use density operator p A into Bob's: p B = 
Na-b{p A )- The two-user quantum broadcast channel 
Na-bc is a quantum channel from sender Alice (A) to 
two independent receivers, Bob (B) and Charlie (C). 
The quantum channel from Alice to Bob is obtained 
by tracing out C from the channel map, i.e., Na-b = 
Trc [Na-bc), with a similar definition for Na-c- We say 
that a broadcast channel Ma-bc is degraded if there ex- 
ists a degrading channel, Af B ^ c , from B to C satisfying 

Ma-C — Nb^c Na-b ■ The degraded broadcast channel 
describes a physical scenario in which for each successive 
n uses of Na-bc Alice communicates a randomly gen- 
erated classical message (m,k) £ (Wb,Wc) to Bob and 
Charlie, where the message sets Wb and Wc have cardi- 
nalities 2 nRB and 2 nRc respectively. The messages (to, k) 
are assumed to be uniformly distributed over (Wb,Wc)- 
Because of the degraded nature of the channel, Bob re- 
ceives both to and k, whereas Charlie only receives k. 

To convey the message (to, k), Alice prepares an n- 
channel-use input state, with density operator k , from 
A n , the tensor product space of her single-use input- 
state alphabet. After transmission through the chan- 
nel, this state results in the bipartite density operator 
P™T = ^ZldP^k) for Bob and Charlie. The re- 
duced density operators for Bob and Charlie, p^ k and 
Pm h respectively, can be found by tracing out the other 
receiver. A (2 nRB ,2 nRc ,n, e) code for this channel con- 
sists of an encoder 

(m,k):(W B ,W c )^A n , (9) 

a positive operator- valued measure (POVM) {A m k} on 
B n and a POVM {A' fc } on C n that satisfy 

Tr(p^ c "(A mfc ®A^)) >l-e, (10) 



for all (m, k) e (Wb,Wc)- Its error probability therefore 
obeys Pi"^ < e. A rate-pair (Rg,Rc) is achievable if 
there exists a sequence of (2 7iRb , 2 nRc , n, e n ) codes with 
e„ — > 0, so that P e ^ for such a sequence. The clas- 
sical capacity region of the degraded quantum broadcast 
channel is then the convex hull of the closure of all achiev- 
able rate pairs (Rb, Rc)- 

The classical capacity region of the two-user degraded 
quantum broadcast channel jVa-bc was recently derived 
by Yard et al. 0], and can be expressed in terms of the 
Holevo information Um. 
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where {pj} is a probability distribution associated with 
the density operators {aj}, and S(p) = — Tr(/51og/S) is 
the von Neumann entropy of the quantum state p. Be- 
cause x mav n °t be additive, the rate region (Rb,Rc) 
of the degraded broadcast channel must be computed by 
maximizing over multiple channel uses. Thus, for n chan- 
nel uses we can achieve the rate region specified by 
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Rc < x(pi,J2P^®c(pf)]/ 7 

S\J2PiPi\iPj 



(12) 
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(13) 



where j = {in, k) is a collective index. The prob- 
abilities {pi} form a distribution over an auxiliary 
classical alphabet T, of size |T|, satisfying \T\ < 
min {\A n \, \B n \ 2 + \C n \ 2 + 1}. The ultimate rate-region 
is computed by maximizing the region specified by 
Eqs. (fT2|) and fl3]), over {pi}, {pj\i}, {pf}, and n, 
subject to the cardinality constraint on \T\. Figure Q] 
illustrates the setup of the two-user degraded quantum 
channel. 



FIG. 1: Schematic of the degraded quantum broadcast chan- 
nel. The transmitter Alice (A) encodes her messages to Bob 
(B) and Charlie (C) in a classical index j by preparing B and 
C in the bipartite state pf ° . The dashed line from B to 
C denotes the existence of a degrading channel, J\fg B _Q, whose 
n-fold tensor product will transform pf into pf for all j. 



IV. NOISELESS BOSONIC BROADCAST 
CHANNEL 

The two-user noiseless bosonic broadcast channel 
Na-bc consists of a collection of spatial and tempo- 
ral bosonic modes at the transmitter (Alice) that inter- 
act with a minimal-quantum-noise environment and split 
into two sets of spatio-temporal modes en route to two in- 
dependent receivers (Bob and Charlie). The multi-mode 
two-user bosonic broadcast channel Na-bc is given by 
® s jVa„-b s c s , where Ma 3 -b 3 C 3 is the broadcast-channel 
map for the sth mode, which can be obtained from the 
Heisenberg evolutions 



b s = y/rk& s + ^/l - rj s e s , 
c s = \/l - r) s a s - ^fq~ s e s , 



(14) 
(15) 



where {a s } are Alice's modal annihilation operators, and 
{bs}, {c e } are the corresponding modal annihilation oper- 
ators for Bob and Charlie, respectively The modal trans- 
missivities {n s } satisfy < rj s < 1, and the environment 
modes {i s } are in their vacuum states. We will limit 
our treatment here to the single-mode bosonic broadcast 
channel, as the capacity of the multi-mode channel can 
in principle be obtained by summing up capacities of all 
spatio-temporal modes and maximizing the sum capac- 
ity region subject to an overall input-power budget us- 
ing Lagrange multipliers, cf. Q, where this was done for 
the capacity of the multi-mode single-user lossy bosonic 
channel. 

The principal result we have for the single-mode de- 
graded bosonic broadcast channel depends on a minimum 
output entropy conjecture (strong conjecture 2, see Ap- 
pendix A). Assuming this conjecture to be true, we will 
show that the ultimate capacity region of the single- mode 
noiseless bosonic broadcast channel (see Fig. ^ with a 
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(ft = |O){0|) 

e 




6 = s/rja + y/l — rje 
1 



c = y/l-rja - y/rje 

Charlie : C {p c } 



FIG. 2: (Color online) A single-mode noiseless bosonic broad- 
cast channel can be envisioned as a beam splitter with trans- 
missivity 77. With 77 > 1/2, the bosonic broadcast channel is a 
degraded quantum broadcast channel, where Bob (B) is the 
less-noisy receiver and Charlie (C) is the more-noisy receiver. 



mean input photon-number constraint 



< N is 



Rb < 9{r)pN), (16) 
Rc < g((l-v)N)-g((l-v)PN). (17) 

Here, < (3 < 1 is a parameter that represents the frac- 
tion of Alice's average photon number that is used to 
convey information to Bob, with remainder to be used to 
communicate information to Charlie. The boundary of 
the broadcast channel's capacity region is traced out by 
varying f3 from to 1. 

It is worth noting, at this point, that our assumption of 
a lossless beam splitter — in Eqs. (fl4|) (|15l) . and Fig. [2] — 
is not essential. In particular, if the coupling coefficients 
from a to b and a to c in Fig. [2] were r\\, and 77 c , respec- 
tively, with < ?7 C < rjb, and rji, + r/ c < 1, then we still 
have a degraded quantum broadcast channel, and, as- 
suming strong conjecture 2 is correct, its capacity region 
is given by Eqs. (| 16[) and (fTT)) , with rjb and r\ c replacing 77 
and 1 — 77, respectively. For simplicity, in all that follows, 
we shall presume that the lossless beam splitter model 
applies. 

The rate region from Eqs. (TT6"|) and (|17[) is additive and 
achievable with single channel use coherent-state encod- 
ing using the distributions 



Pr(t) 



PA\r{a I t) 



1 



exp 



w 

N 



exp 



Nf3 



(18) 



(19) 



The first step toward proving that Eqs. (fT6|) and (fl~7|) 
do indeed specify the bosonic broadcast channel's capac- 
ity region is to show that Eqs. (JTSj) and (fl9|) achieve these 
rates. It is straightforward to show that if 77 > 1/2, the 
bosonic broadcast channel is a degraded quantum broad- 
cast channel, in which Bob's is the less-noisy receiver 



and Charlie's is the more- noisy receiver. To do so we 
merely recognize that, when 77 > 1/2, Charlie's reduced 
density operator can be obtained from Bob's by applying 
{i>i ■ 1 < i < 71} to the input of a lossless beam splitter 
that has transmissivity 77' = (1 — T])/r] to output modes 
{di : 1 < i < n}, and whose other input port is driven 
by vacuum-state modes. The Yard et al. capacity region 
in Eqs. (|12|) and (fl~3|) requires finite-dimensional Hilbcrt 
spaces for the channel's input and outputs. Nevertheless, 
we will use their result for the bosonic broadcast chan- 
nel, which has an infinite-dimensional state space, by ex- 
tending it to infinite-dimensional state spaces through a 
limiting argument [151 ]. The n = 1 rate-region for the 
bosonic broadcast channel using a coherent-state encod- 
ing is thus: 

Rb < J p T {t)S^J PAlT {a\t)\^ja)(^ja\da^jdt (20) 
Rc < S^J p T {t)p A \T(a I t)\y/l -va)(y/l =rja\ dad?j 



(21) 



Pr(t)Syj p A \r{a \ t) 
x I -y/l — rja)(\/l — r/a\ daj dt, 



where we need to maximize the bounds for Rb and Rc 
over all joint distributions Pt{^)Pa\t{P- I t) subject to 
(|a| 2 ) < N. Note that A and T are complex-valued ran- 
dom variables, and the second term in the Rb bound 
(|12[) vanishes, because the von Neumann entropy of a 
pure state is zero. Substituting Eqs. (fl"8|) and (fT9f into 
Eqs. (|20|) and ((2"l"|) shows that the rate region in Eqs. (fl6|) 
and (fT7|) is achievable with single-use coherent-state en- 
coding. 

For the converse, assume that the rate pair (Rb,Rg) 
is achievable. Let {(m, k)} and the POVMs {A m k}, 
{A' k } comprise any (2 nRB , 2 nRc , n, e„) code in the achiev- 
ing sequence. Suppose that Bob and Charlie store 
their decoded messages in the classical registers Wb 
and Wc respectively. Let us use pw B ,w c { m ^) — 
Vw B (m)p Wc ( k ) = 2- nHB 2- ,lRc to denote the joint prob- 
ability mass function of the independent message regis- 
ters Wb and Wc- As (Rb,Rc) is an achievable rate-pair, 
there must exist e' n — > 0, such that 



nR c = H{W C ) 

< I{W c ;W c )+ne' n 

< x{Pw c ( k )-,Pk n ) "+ 



ne r 



(22) 



where I(W C ;W C ) = H{W C ) - H(W C \ W c ) gives the 
Shannon mutual information in terms of the Shannon 
entropy H = —^2 k Pk^-OE,{Pk) for a probability distribu- 
tion {p k }, and p^" = J2 m Pw B ( m )pZk- Thc second line 
follows from Fano's inequality [l6[ and the third line fol- 
lows from Holevo's bound [17|- Similarly, for e" — > 0, we 
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can bound uRb as follows: 

nR B = H(W B ) 

< I(W B ;W B ) + ne'^ 

< ^Pw c (k)x(Pw B (m),p? n " k ) +ne' r l, (23) 

fc 

where the three lines above follow from Fano's inequality, 
Holevo's bound and the concavity of Holevo information. 

To complete the converse proof, we need only show 
that there exists a < (3 < 1, such that 

^2pwc(k)x(Pw B (m),p? n " k ) < ng(r){3N), (24) 
fc 

X(pwc(k),pf) < ng((l - V )N) - ng((l - V )pN). (25) 

From the non-negativity of the von Neumann entropy 
S {pZ,k)i it; follows that J2kPw c ( k )x(Pw B (m), p^ k ) < 

J2kPw c m Pw B (m)p^ k ), as the second term of 

the Holevo information above is non-negative. Because 
von Neumann entropy is subadditive, and the maximum 
von Neumann entropy of a single-mode bosonic state 
with (fit a) < N is given by g(N), we have that 

n 

< S[pf) < Y,9(vN ke ) < ng(riN k ) , (26) 

where, N k = Yle=i Nk e /n, and N kl is the mean pho- 
ton number of the £th channel use for the state p k " = 

J2 m Pw B ( m )pm.k- Therefore, because g(0) = and g(x) 
is monotonically increasing for x > 0, we see that for 
each k G Wc there is a < [3 k < 1 such that 

s(pf)=ng{r ] p k N k ). (27) 

We know that N is Alice's maximum average photon 
number per channel use, where the averaging is over the 
entire codebook. Thus, the mean photon number of the 
n-use average codeword at Bob, p B " = J2kPw c (k)Pk > 
is rjN . Hence, we get 

< Y,PWc{k)s(pT) < S( P B ") < ng ( V N) , (28) 
fc 

where the second inequality follows from the convexity of 
von Neumann entropy. Again invoking the monotonicity 
of g(x) we have that there is a < j3 < 1, such that 
Y.kPw c {k)S(pf l ) = ng(r](3N) whence 

^2pw c ( k )x(Pw B (m),p^ k ) <ng{rif3N). (29) 

fc 

This proves the first inequality that we need for the ca- 
pacity region's converse statement. 



To prove the second inequality needed for that con- 
verse, we start from Eq. (|2T|) and use strong conjecture 2 
(see Appendix A) to get 

s(f>Z ) >ng({l-rj)(3 k N k ). (30) 

Next, we use the uniform distribution Pw c ( k ) — 2~ nRc 
to obtain 

2~ nRc 9 (vPkNk) = g (v(3N) ■ (31) 

fc 

Using (f3"Tj) . the convexity of g(x) and r\ > 1/2, we have 
shown (see Appendix B) that 

]T 2- nR °g ((1 - V )faN k ) > 9 {(l- r,)W) ■ (32) 
fc 

From Eq. Q32p . and Eq. pop summed over k, we then 
obtain 

Y,Pw c (k)s[pf ) >ng((l- V )(3N). (33) 
fc 

Finally, writing Charlie's Holevo information as 

x(Pw c { k ),Pk' 1 ) = 

V & / fc 

< ng{{l-ri)N)-Y,PWo{k)s(pf ), (34) 
fc 

we can use Eq. (|3"3"|) to get 

x(p Wc (fc),/3r)<n 5 ((l-T?)iV)-n 5 ((l-77)/3iV), (35) 

which completes proof of the converse, given that strong 
conjecture 2 is true. 

V. DISCUSSION 

Without a proof of strong conjecture 2, we cannot as- 
sert that Eqs. (fTrJ)) and (fT7| define the capacity region 
of the bosonic broadcast channel. However, because the 
rate region specified by these equations is achievable — 
with single-use coherent-state encoding — we know that 
they comprise an inner bound on that capacity region. 
In this regard it is instructive to examine how the rate 
region defined by Eqs. (| 16[) and (| 1 7[) compares with what 
can be realized by conventional, coherent-detection op- 
tical communications. Suppose Alice sends a coherent 
state, to the beam splitter shown in Fig. [5] Bob 
and Charlie will then receive coherent states \yfrja) and 
— rj a), respectively. Moreover, if Bob and Char- 
lie employ homodyne-detection receivers, with local os- 
cillator phases set to observe the real-part quadrature, 
their post-measurement data will be y/rjHe(a) + vb for 



6 



Bob and \J\ — 77Re(a) + vc for Charlie, where vb and 
vc are independent, identically distributed, real- valued 
Gaussian random variables that are zero mean and have 
variance 1/4 [l8[ . Similarly, if Bob and Charlie employ 
heterodyne-detection receivers, their post-measurement 
data will be y/rja + zb and yjl — rja + zc, where zb and 
Zc are independent, identically-distributed, complex- 
valued Gaussian random variables that are zero mean, 
isotropic, and have variance 1/2 d8j. These results im- 
ply that the r\ > 1/2 bosonic broadcast channel with 
coherent-state encoding and homodyne detection is a 
classical degraded scalar-Gaussian broadcast channel, 
whose capacity region is known to be fl9j 



Rb < ^ log(l + 



Rc < 2 lo g( 1 



4(1 - r,)(l - 0)N) \ 
1 + 4(1 - rj)PN J ' 



(36) 
(37) 



for < j3 < 1. Likewise, the r\ > 1/2 bosonic broadcast 
channel with coherent-state encoding and heterodyne de- 
tection is a classical degraded vector- Gaussian broadcast 
channel, whose capacity region is known to be 

Rb < log(l + TiPN) (38) 

for < < 1. In Fig. [3] we compare the capacity re- 
gions attained by a coherent-state input alphabet using 
homodyne detection, heterodyne detection, and optimum 
reception. As is known for single-user bosonic communi- 
cations, homodyne detection performs better than het- 
erodyne detection when the transmitters are starved for 
photons, because it has lower noise. Conversely, hetero- 
dyne detection outperforms homodyne detection when 
the transmitters are photon rich, because it has a factor- 
of-two bandwidth advantage. To bridge the gap between 
the coherent-detection capacity regions and the ultimate 
capacity region, one must use joint detection over long 
codewords. Future investigation will be needed to de- 
velop receivers that can approach the ultimate commu- 
nication rates over the bosonic broadcast channel. 

Recently, Vishwanath et al. [1] established the du- 
ality between the dirty-paper achievable rate re gion — 
recently proved to be the ultimate capacity region [251] — 
for the classical multiple-input, multiple-output (MIMO) 
Gaussian broadcast channel and the capacity region of 
the classical MIMO Gaussian multiple-access channel 
(MAC). Their duality result states that if we evaluate 
the capacity regions of the MIMO Gaussian MACs — with 
fixed total received power P and channel-gain values 
over all possible power allocations between the users, 
the corners of those capacity regions trace out the ca- 
pacity region of the MIMO Gaussian broadcast channel 
with transmitter power P and the same channel-gain val- 
ues. Of course, the bosonic broadcast channel and the 
bosonic multiple-access channel satisfy this duality when 




FIG. 3: (Color online) Comparison of bosonic broadcast chan- 
nel capacity regions, in bits per channel use, achieved by 
coherent-state encoding with homodyne detection (the capac- 
ity region lies inside the boundary marked by circles), hetero- 
dyne detection (the capacity region lies inside the boundary 
marked by dashes), and optimum reception (the capacity re- 
gion lies inside the region bounded by the solid curve), for 
n = 0.8, and N = 1, 5, and 15. 



they employ coherent-state encoding and coherent de- 
tection, because under these conditions these quantum 
channels reduce to classical additive Gaussian-noise chan- 
nels. However, it turns out that the capacity region of 
the bosonic broadcast channel using coherent-state in- 
puts and optimum reception is not equal to of the en- 
velope of the MAC capacity regions using coherent-state 
inputs. The capacity region of the bosonic MAC using 
coherent-state inputs was first computed by Yen 0]. In 
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Rb 




FIG. 4: (Color online) Comparison of bosonic broadcast and 
multiple-access channel capacity regions, in bits per channel 
use, for n = 0.8, and N = 15 with coherent-state encoding. 
The dashed curve is the outer boundary of the broadcast ca- 
pacity region. It lies below the solid curve, which is the outer 
envelope of the MAC capacity regions. 

Fig. Q] we compare the envelope of coherent-state MAC 
capacities to the capacity region of the coherent-state 
broadcast channel. This figure shows that with a fixed 
beam splitter and identical average photon number bud- 
gets, more collective classical information can be sent 
when the beam splitter is used as a multiple-access chan- 
nel as opposed to when it is used as a broadcast channel 
if coherent-state encoding is employed. 



into the channel output. The study of minimum output 
entropy yields important information about channel ca- 
pacities, viz., an upper bound on the classical capacity 
derives from a lower bound on the output entropy of mul- 
tiple channel uses, see, e.g., [5j. Also, the additivity of 
the minimum output entropy implies the additivity of the 
classical capacity and of the entanglement of formation 

In this appendix we first briefly review previous work 
on a minimum output entropy conjecture that arose in 
conjunction with the channel capacity analysis of the 
single-user bosonic channel with additive Gaussian noise 

H,ll> US-SI- Thcn wc wiU turn our attention to the 
minimum output entropy conjecture that we employed 
in our capacity theorem for the degraded bosonic broad- 
cast channel. Both conjectures have weak (single-use) 
and strong (multiple- use) versions. 

Let a and b denote the two input modes of a lossless 
beam splitter of transmissivity 77, that has output modes 
c = ^/rj a + y/1 — i] b and d = ^/l — rj a — Jr\ b. In [6| , the 
following minimum output entropy conjecture was made. 

Conjecture 1 — Let the input mode b be in a thermal 
state with average photon number K (hence von Neu- 
mann entropy g{K) ). Then the von Neumann entropy of 
the output mode c = ^/rja + y/1 — nb is minimized when 
the input mode a is in the vacuum state. The resulting 
minimum von Neumann output entropy is g((l — T])K). 

The above conjecture is a special case of the following 
strong (multi-mode) version whose proof would establish 
the ultimate capacity of the single-user bosonic channel 
with thermal noise. 
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APPENDIX A: MINIMUM OUTPUT ENTROPY 
CONJECTURES 

In general, the evolution of a quantum state resulting 
from the state's propagation through a quantum com- 
munication channel is not unitary, so that a pure state 
loses some coherence in its transit through the channel. 
The minimum von Neumann entropy at the channel's 
output provides a useful measure of the channel's ability 
to preserve the coherence of its input state. In partic- 
ular, the output entropy associated with a pure state 
measures the entanglement that such a state establishes 
with the environment during propagation through the 
channel. Because the state of the environment is not ac- 
cessible, this entanglement is responsible for the loss of 
quantum coherence and hence for the injection of noise 



Strong Conjecture 1 Let the input modes 

b = [bi 62 • ■ • b n } T be in a product state of n ther- 
mal states with total von Neumann entropy ng(K). 
Then the von Neumann entropy of the output modes 
c = [ ci C2 • ■ • c n ) T = ^/rja + s/1 — r\ b is minimized 
when the input modes a = [ &i S2 • • • a n ] T are in their 
vacuum states. The resulting minimum von Neumann 
output entropy is ng((l — n)K). 

Neither strong conjecture 1 nor its weak (single-use) 
form have been proven yet, but considerable evidence in 
support of their validity has been developed. For exam- 
ple, strong conjecture 1 has been shown to be true when 
the input states are restricted to be Gaussian (23J. It 
has also been proven that the vacuum state provides a 
local minimum for the output entropy [6(. Strong con- 
jecture 1 has been shown to be true when Renyi entropy 
of integer order > 2 is employed in lieu of von Neumann 
entropy [2l| . Similarly, conjecture 1 has been proven 
when Wehrl entropy — the continuous Boltzmann-Gibbs 
entropy of the Husimi probability function [24j — is used 
instead of von Neumann entropy [2l[. Additional evi- 
dence in support of conjecture 1 can be found in Q. 



8 



In proving the converse to the bosonic broadcast 
channel's capacity theorem we assumed the validity of 
the following conjecture. 

Strong Conjecture 2 — Let the input modes a be in 
their vacuum states, and let the von Neumann entropy 
of the input modes b be ng(K). Then, putting the b 
modes in a product state of mean-photon-number K 
thermal states minimizes the von Neumann entropy of 
the output modes c = ^/rj a + y/1 — 77 6. The resulting 
minimum von Neumann output entropy is ng((l — rj)K). 

The weaker, single-use version of this conjecture is 
also of interest. 

Conjecture 2 — Let the input mode a be in its vacuum 
state, and the let the von Neumann entropy of the input 
mode b be g{K). Then the von Neumann entropy of the 
output mode c — ^/rj a + ^1 — rj b is minimized when b 
is in a thermal state with average photon number K. 
The resulting minimum von Neumann output entropy is 
<?((!- 

We have yet to develop proofs for either strong con- 
jecture 2 or conjecture 2. In the rest of this appendix 
we will present evidence that supports their validity. To- 
ward that end, we first show that strong conjecture 2 is 
true when Wehrl entropy is used instead of von Neumann 
entropy. 



1. Strong Conjecture 2 for Wehrl Entropy 

Strong Conjecture 2: Wehrl — Let the input modes 
a be in their vacuum states, and let the Wehrl entropy of 
the input modes b be n(l + hi(K + 1)). Then, putting the 
b modes in a product state of mean-photon-number K 
thermal states minimizes the Wehrl entropy of the output 
modes c = y/rja + \J1 — rj b. The resulting minimum 
Wehrl output entropy is n(l + In (K(l — rf) + 1)). 

Proof — The Wehrl entropy for an n-mode density op- 
erator p is 



namely 



W(p) 



Qp( M ) In KQp( M )]d 2 >, 



(Al) 



where Qp(fi) = (/j|/5|/x)/7r n , with \fi) a coherent state, 
is the Husimi distribution, i.e., the joint probability den- 
sity function for multi-mode heterodyne detection. The 
Wehrl entropy provides a measurement of the state p in 
phase space and its minimum value is achieved on coher- 
ent states [24j . 

Our proof of strong conjecture 2 for Wehrl entropy re- 
lies on the antinormally-ordered characteristic function, 
Xa"(Oj associated with the n-mode density operator p, 



X P A 



•(0=tr(p a e-< t& e< T&t ), (A2) 

where £ = [ Ci C2 • • • Cn ) T is a column vector of 
complex numbers, = [ Ci C2 '"" Cn L an d a* = 
[d\ a\ ■ ■ ■ ajj ] T . The antinormally-ordered character- 
istic function and the Husimi function are a 2n-D Fourier 
transform pair: 



X P A (0 = / Qp(M)e 



(A3) 



M = -57 / xtoe-^^C (A4) 



As the two n-use input modes a and b are in a product 
state, Eq. (|A2|I implies that the output-state characteris- 
tic function is a product of the input-state characteristic 
functions with appropriately scaled arguments, 



x p a (0 = x p a (VvC)x P A b (V^vO- (A5) 

From Eq. (|A5[) . the multiplication-convolution and scal- 
ing properties of Fourier-transforms pairs, and the fact 
that a is in the n-mode vacuum state, we find that 



1 



(7rn) n 



1 



(A6) 



where * denotes convolution. 

Suppose that the state of the input modes b is a prod- 
uct of thermal states, each with mean photon number K, 
i.e., 



Pb 



(-L f e-^ 2 / K \a)(a\d 2 a 



\ttK 



(A7) 



The Wehrl entropy for the b modes is then 

W(p b ) = n(l + hx(K + 1)), (A8) 

which satisfies the hypothesis of strong conjecture 2 for 
Wehrl entropy. Using Eq. (|A6[) . we can now show that 
the Husimi function and the Wehrl entropy for the state 
of the output modes c are 



-M a /{i+0- v )K) 



(A9) 



Wl + Cl-r?)^))"' 
W( Pc ) = n(l + ln(^(l-7 ? ) + l)), (A10) 

providing an upper-bound to the minimum Wehrl output 
entropy. 

To show that the expression in Eq. (|A8|) is also a lower 
bound for the Wehrl output entropy, we use Theorem 6 
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of [26| , which states that for two probability distributions 
f(fJt) and h([i) over n-dimensional complex vectors fj, we 
have 

W((f * h)(ti)) > XW(f(ti)) + (1 - X)W(h(fi)) 

-nAlnA-n(l- A)ln(l-A), (All) 

for all < A < 1, where the Wehrl entropy of a prob- 
ability distribution is found from Eq. (|A1|) by replacing 
Qp(fi) with the given probability distribution. Choosing 



/(M) 



i 



(A12) 
(A13) 



(l-77)"^ pb 

we get 

W(p c ) > nA(l + ln?7) 

- nAlnA - n(l — A)ln(l - A). (A14) 

The Wehrl entropy of a scaled distribution 
(l/x) n Q(fi/ \fx) is easily shown to satisfy 



W 



1 



Q 



W(Q(iJ,))+nlnx, 



(A15) 



for any x > 0. From Eqs. (|A15|) and (|A14[) we then 
obtain 

W(p c ) > nA(l + In 77) + (1 - A) (W(p b ) 
+ nln(l - ij)) - nAln A - n(l - A) ln(l — A) 
= nA(l + In 77) + n[l - A) (1 + ln(K + 1) 
+ ln(l-ry)) - nAlnA - n(l - A)ln(l - A) 
= n(l + ln(if(l-77) + l)). (A16) 

The last equality used \ = i] / (i] + (K + — 77)), which 
satisfies < A < 1 for all ij and K. Therefore the mini- 
mum Wehrl entropy of the output modes c has the lower 
bound n(l + \n(K(l — rf) + 1)). Because this lower bound 
coincides with the upper bound, derived earlier, we know 
that it is indeed the minimum Wehrl output entropy, and, 
moreover, that this minimum is achieved by a product 
thermal-state p b with mean photon number K per mode. 



2. Strong Conjecture 2 for Gaussian-State Inputs 



In this section we prove that stron 
2 are equivalent when all inputs are 
Gaussian states. Because strong con 
proven for Gaussian-state inputs [23| 
implies the truth of strong conjecture 

With no loss of generality we shall 
tion to zero-mean Gaussian states. A 



; conjectures 1 and 
restricted to be in 
iecture 1 has been 
|, this equivalence 
2 for such inputs, 
restrict our atten- 
zero-mean n-mode 



Gaussian state is completely characterized by its 2n x In 
correlation matrix 



Ro = 



it 



[(*) 



T iTl \ _ 



(aa 



(aa T ) 



(A17) 

where I„ is the n x n identity matrix and * denotes 
component-wise complex conjugation. Of particular im- 
portance, for what will follow, is the symplcctic diago- 
nalization of R a and the consequences thereof. 

Let the n modes represented by a be in a zero-mean 
Gaussian state with correlation matrix R a . We will show 
that there exists a 2n x 2n complex-valued symplectic 
matrix S such that 



R a = SAS f , 

where is the conjugate transpose of S, 
StQS = SQS t = Q, 

and 



A = diag[ Ai + 1 • 
Equation (|A19|) . with 

Q 



A„ + l Ai 



I» 
-I„ 



A n 



(A18) 



(A19) 



(A20) 



(A21) 



is the condition that makes S symplectic. The {A^} 
are the symplectic eigenvalues of R a , which are all non- 
negative because R a is positive semidefmite. 

To establish the preceding symplectic diagonalization 
of R a , we use Williamson's symplectic decomposition 
theorem on the symmetrized (real-valued) correlation 
matrix for the quadratures, di = Re(d) and 0,2 = Im(d) 
[13]. Equations (|A18|) - (|A20|) are then obtained by con- 
verting this quadrature correlation-matrix decomposition 
into the annihilation operator correlation matrix via the 
linear transformation 



U 



1 n zl n 



(A22) 



The value of the symplectic decomposition lies in es- 
tablishing a linear relationship between the modes d, 
which are in a given zero-mean Gaussian state, to a 
new set of modes that are in independent thermal states 
whose average photon numbers are given by the sym- 
plectic eigenvalues. In particular, for d in an arbitrary 
n-mode zero-mean Gaussian state with correlation ma- 
trix R a , we have that 



it 



(A23) 



accomplishes this transformation, where S 1 = QS^Q. 
The n modes represented by d are in a zero-mean Gaus- 
sian state with a correlation matrix that is easily found 
to be 



Rrf 



(A24) 
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Thus, di has average photon number (d\di) = Aj, for 1 < 
i < n. Furthermore, the di modes are all uncorrelated — 
because A is diagonal — so that each mode can be rep- 
resented as an isotropic Gaussian mixture of coherent 
states, and the joint state is the tensor product of n such 
states. 

The symplectic transformation in (|A18|) is canonical, 
i.e., it preserves the commutation relations. Thus it can 
be implemented with a unitary operator U, satisfying 
UU^ = &U = I dHH. From tnis w e see that the von 
Neumann entropy of the Gaussian-state a modes is iden- 
tical to that of the thermal-state d modes from (|A23|) . 
We are now ready to address the central concern of this 
section, namely showing that strong conjecture 2 is true 
when the input states are restricted to be Gaussian. 

Theorem 1. Strong conjecture 1 and strong conjecture 2 
are equivalent when the input fields are restricted to be in 
Gaussian states. 

Proof. Consider the the vector input-output relation 



= ^l 



a 



(A25) 



with the a and b modes being in independent quantum 
states that are zero-mean Gaussians. 

First let us use the truth of strong conjecture 1 to 
show that strong conjecture 2 is also true. Under the 
premise of strong conjecture 2, we take the a modes to 
be in their vacuum states, and the b modes to be in 
a zero-mean Gaussian state with correlation matrix R5 
and von Neumann entropy ng{K). (No loss in general- 
ity ensues from the restriction that the b modes be in 
a zero-mean state, because von Neumann entropy is in- 
variant to state displacement.) Because the inputs are in 
Gaussian states, minimizing the von Neumann entropy 
of the output modes c reduces to finding the correlation 
matrix R& that minimizes this entropy. Let R& = SAS^ 
be the symplectic diagonalization of R&. We can then ex- 
press the input modes b as a symplectic transformation 
on another set of n-modes, d, 



(A26) 



whose correlation matrix is R^ = A. Furthermore, we 
have that 

n n 

S(p d ) =5>03 d J = 5>(A,) = S(p b ) = ng(K). 

i=l i=l 

(A27) 



b 




d 




= S 




b\ 




d\ 



Substituting Eq. (|A26|) into (|A25)) and using some lin- 
ear algebra, we get 



s v^s- 1 



a' 



d 



(A28) 




c=Jrjd+Jl-r] b 



S (S f QS - Q) 



T 



^ Product state, zero-mean, 
■ 4 Gaussian, thermal: 



,/,/) A,\.£</(Aj ng(K) 




re-mode vacuum 
state 



Product state, zero-mean, 
Gaussian, thermal: 



FIG. 5: (Color online) Schematic of the beam splitter chan- 
nel with zero-mean Gaussian- state inputs, and the equivalent 
channel after transformations that preserve von Neumann en- 
tropies. 



A schematic corresponding to this equation is shown in 
the bottom panel of Fig. In particular, the beam 
splitter channel governed by Eq. (|A25[) and the Gaussian 
states we have assumed for a and b are equivalent to what 
we have shown in the top panel of Fig. [51 We know that 
symplectic transformations do not affect von Neumann 
entropy. Thus minimizing the von Neumann entropy of 
the c modes by choice of the correlation matrix R& in the 
top panel of Fig. [5] is equivalent to minimizing this out- 
put entropy by choice of the 2n x 2n symplectic matrix 
S and the symplectic eigenvalues {Ai>0:l<i<n}, 
subject to the constraint that X^i=i3(^«) = n ff(^0? m 
the lower panel of that figure. 

Suppose that we have a set of symplectic eigenvalues 
that satisfy the constraint. Then, via strong conjecture 1, 
the von Neumann entropy the c modes is minimized when 
the a modes in the lower panel of Fig. [5] are in their vac- 
uum states. However, because the a modes are already in 
this state, an optimizing symplectic transformation S" 1 
is the identity matrix This result is independent of 
the particular values of the {A^}, but the entropy of the c 
modes still depends on our choice of symplectic eigenval- 
ues. In particular, when the a modes are in their vacuum 
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states and S 1 = l2 n , the von Neumann entropy of the 
c modes is 



(A29) 



To minimize this output entropy, by choice of the {Aj} 
we employ a Lagrange multiplier approach: 



min > o(vi — ri A;) 

{A»>0:E? =1 s(A i )=n ff (i<:)} ^ 



(A30) 



Differentiating Eq. (|A30[) with respect to the {A^} and £ 
yields 



£ = 77~\ \ , for 1 < « < n (A31) 

.9 (Ai) 

n 

X>(Ai) = n0(tf) , (A32) 

i=l 

which implies that choosing Ai = if, for 1 < i < n, 
minimizes the output entropy subject to the constraint 
[3fj| . The minimum output entropy is then ng(y/l — rjK), 
and it is achieved when the n-mode Gaussian input state 
is an n-mode thermal product state with (b\bi) = K for 
1 < i < n. This completes the demonstration that strong 
conjecture 1 implies strong conjecture 2 for Gaussian- 
state inputs, and because strong conjecture 1 is known 
to be true for Gaussian-state inputs, we have proven that 
strong conjecture 2 is also true for such inputs. To com- 
plete the proof of Theorem 1, we must still show that 
strong conjecture 2 implies strong conjecture 1 when the 
input states are Gaussian. 

Assume that strong conjecture 2 is true, and let the in- 
put modes b be in a product state of n zero- mean thermal 
states each with von Neumann entropy g{K), as shown 
in the top panel of Fig. [6] With no loss of generality 
we can take the input modes a to be in a zero-mean 
pure Gaussian state, i.e., a is in an n-mode vacuum or 
squeezed- vacuum state. Performing the symplectic diag- 
onalization R a = SAS\ we write 



(A33) 



where = A. Because this transformation preserves 
von Neumann entropy, we know that must be a zero- 
mean pure Gaussian state with no phase-sensitive corre- 
lation. The only such state is the n-mode vacuum state. 
We can then perform similar algebraic manipulations to 
the beam splitter relation in Eq. (|A25[) to get, 



a 




d 




= S 






d\ 



b 



(A34) 



n-mode vacuum ^ 

state <4» 



Product state, zero-mean, 
Gaussian, thermal: 

#6A = 2jfft 




c=Jrja+Jl-r] b 



Product state, zero-mean, 
Gaussian, thermal: 



i 3 



Kb, 



FIG. 6: (Color online) Schematic of the beam splitter chan- 
nel with zero-mean Gaussian- state inputs, and the equivalent 
channel after transformations that preserve von Neumann en- 
tropy. 



as shown in the lower panel in Figure [6l 

Minimizing the von Neumann entropy after the sym- 
plectic transformation at the output port in the lower 
panel of Fig. [5] is equivalent to minimizing the entropy 
before that transformation. Thus our objective is to de- 
termine the 2n x 2n symplectic matrix S _1 that mini- 
mizes the von Neumann entropy before the output-port 
symplectic transformation. Because the a modes are in 
their vacuum states and the modes applied to the beam 
splitter's other input port have von Neumann entropy 
ng(K), strong conjecture 2 tells us that the latter in- 
put should be in an n-mode thermal product state, with 
average photon number K per mode, to achieve the min- 
imum output entropy. But b is already in this state, 
so an optimizing symplectic transformation is therefore 
the identity, S _1 = Ian.. This allows us to conclude that 
putting the a modes in their vacuum states minimizes 
the entropy of the c modes when the b modes are in an 
n-mode product of thermal states each with average pho- 
ton number K, thus demonstrating that strong conjec- 
ture 2 implies strong conjecture 1 when the input states 
are Gaussian. □ 
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FIG. 7: (Color online) Plots of the 1-to-l function g(x) and 
its inverse h(y), showing y = h^ 1 [rjh(y)) = g { r gg~ 1 {y)) and 

Kv) = y~y'- 



APPENDIX B: A PROPERTY OF g(x) 

For the converse proof given in Sec. IV, we need to 
show that for non-negative real numbers {xk : 1 < k < 
n}, and < r\ < 1, if Xq is defined by 



E 



g(xk) 



k=l 



then 



E 

k=l 



(Bl) 



(B2) 



where g(x) = (1 + x) log(l + x) — xlog(x). 

Because g(x) is a 1-to-l function, it has an inverse 
function h(y) = g~ 1 (y), such that if y = g{x) then x = 
h(y). Let y/j = g(xk), for 1 < k < n. For every y > 0, 
define y' = (rjh(y)) = g (vg~ 1 (y)) an( i define = 
j/ — j/', as shown in Fig. [7l Using this notation what we 
are trying to prove becomes the following. Given that 



I " 

II ' 



(B3) 



fe=i 



show that 



1 - 



(B4) 



fc=i 



By straightforward differentiation, we can show that 
d 2 l(y)/dy 2 < which implies that 



n 1 n 

yo-y' >-^yk- ~E^> ( B5 ) 
fc=i fe=i 

from the definition of /(j/). Using Eq. (|B3j) we then get 

1 " 

fc=i 

which completes the proof. 
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